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We consider the nonrelativistic interacting Bose gas at zero and finite temperature. Using the 
1/iV-expansion, we derive expressions for the free energy density and the number density to next- 
to-leading order in 1/N. Outside the critical region and at weak coupling, our calculations reduce 
to the well-known mean-field results for the dilute Bose gas. We also rederive the nonperturbative 
critical density for Bose condensation which was first calculated by Baym, Blaizot, Holzmann, Laloe, 
and Vautherin. 
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I. INTRODUCTION 



The dilute Bose gas has a very long history dating back to 
the classic paper by Bogoliubov in 1947 Q]. At zero tempera- 
ture, the quantum loop expansion is essentially an expansion 
in powers of the gas parameter pa 3 , where p is the number 
density and a is the s-wave scattering length. The leading 
quantum corrections to the mean-field results for the energy 
density and the speed of sound were calculated many years 
ago by Lee and Yang 2] . 

Finite-temperature corrections to the pressure of a dilute 
Bose gas were first calculated by Lee and Yang [3J. By per- 
forming an expansion of the pressure about zero temperature, 
they showed that the leading correction term goes as T 4 . This 
reflects the fact that the thermodynamics at low temperature 
is completely dominated by the linear part of the Bogoliubov 
spectrum. 

Most approaches to the thermodynamics of the dilute Bose 
gas arc mean- field approximations SBIi0EH[I3,[H 
fl3 . Il3l IT^jl . Common to all mean-field calculations is that 
they predict a critical temperature T c for Bose condensation 
which is the same as that of an ideal Bose gas |l5j| . More- 
over some of the mean-field approximations such as the Popov 
appoximation predict a first-order phase transition. This is 
in disagreement with universality-class arguments which tell 
us that the NR Bose gas goes through a second-order phase 
transition since the it belongs to the universality class of the 
3-dimensional 0(2)-symmetric spin model. Another prob- 
lem with some of the mean-field approximation such as the 
Hartree-Fock-Bogoliubov (HFB) approximation, is that the 
dispersion relation is gapped 151 [16] . This is in conflict with 
the Hugenholtz-Pines theorem \VA Il8| which states that the 
spontaneous breaking of the U(l) symmetry in the nonrel- 
ativistic (NR) Bose gas gives rise to a gapless mode. Note 
however, that this problem with the HFB approximation has 
recently been solved by Yukalov and Kleinert by introducing 
a separate chemical potential for the condensate 14) (sec also 
Ref. [l9l| for a modfied gapless Hartree-Fock approach in the 
context of relativistic 0(A^)-symmetric theories). 

Examples of calculations beyond mean-field theory are the 
renormalization group calculations of the spin-zero Bose gas 
of Refs. [2^, CHI E3 • Such calculations typically sum up 



classes of diagrams from all orders of perturbation theory and 
they show that the dilute Bose gas undergoes a second-order 
phase transitions in accordance with expectation based on 
universality [U \z&. Other calcuations |2 it |25|| are based 
on the two-particle-irreducible (2PI) diagram effective-action 
approach first developed by Luttinger and Ward 26] and by 
Baym |27| in the context of nonrelativistic fermions. Com- 
bined with the 1/N expansion, the 2PI effective action is par- 
ticular suited to deal with nonequilibrium phenomena |25j|. 
and has also been used successfully to calculate critical ex- 
ponenents [2S||. Very recently, the equation of state of the 
dilute Bose gas was calculated by Pilati et al |29l| using the 
Path Integral Monte Carlo method. 

The problem of calculating the critical temperature for 
a dilute Bose gas has been around for half a century (See. 
Refs. [l6, 30] for a review), butthe issue been was first settled 
6-7 years ago by Baym et al |3lll32jl. The Hartree-Fock (HF) 
approximation which amounts to perturbation theory to lead- 
ing order in the scattering length merely implies a redefinition 
of the chemical potential and hence predicts no shift in the 
critical temperature. Going beyond the HF approximation 
using perturbation theory, one is facing infrared divergences 
which become more severe as one goes to higher orders. Thus 
perturbation theory breaks down close to T c which shows that 
the physics of the phase transition is inherently nonperturba- 
tive. Using effective-field-theory methods, it was shown rigor- 
ously by Baym et al (3ll |3^| that the parametric dependence 
of the change in the critical temperature, AT C , is linear in the 
scattering length a in the dilute limit: 
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where c is a constant, and T° and p2 are the critical temper- 
ature and critical density of the ideal Bose gas, respectively. 
The coefficient c has been determined by analytical as well 
as numerical means in recent years. This includes effective- 
field theory f3l|. I3I [SI . high -precision Monte-Carlo calcula- 
tions |34l l3Bl I36L l37j , simulations of classical field theory in 
the microcanonical ensemble l3g||. 1 /W-expansions |39l |40| , 
variatio nal p erturbation theory |41|. and the linear delta ex- 
pansion Him [3. 

The fact that mean-field theories and perturbation the- 
ory fail in the critical region, warrants a nonperturba- 
tive approach. In the present paper, we apply the 1/N- 
expansion [45I . liril] to the interacting NR Bose gas at finite 
temperature. In contrast with e.g. the HFB approximation, 
this approach satisfies the Hugenholtz-Pines theorem order by 



order in 1/TV. This expansion is a nonperturbative method 
that sums up infrared-divergent diagrams from all orders of 
perturbation theory. For example, a leading-order calculation 
sums up all daisy and superdaisy diagrams. Moreover, in or- 
der to get results beyond mean field, one must typically go 
to nex-to- leading order (NLO) in the 1/TV expansion. Follow- 
ing the approach in Ref. |47|. we derive the thermodynamic 
quantities at finite temperature to next-to-leading order in 
the 1/TV-expansion. 

The paper is organized as follows. In Sec. II, we briefly 
discuss the nonrelativistic Bose gas and the 1/TV-expansion. 
In Sec. Ill, we discuss the thermodynamic potential and the 
gap equations. We also derive the leading-order and next-to- 
leading-order results. We summarize in Sec. IV. 



II. DILUTE BOSE GAS AND EFFECTIVE 
ACTION 

The Euclidean Lagrangian for a nonrelativistic Bose gas 
with N species is 

(2) 

where i = 1,2, N, ipi is a complex field, fi is the chemical 
potential 1 and M is the mass of the bosons. In the following 
we will be using units where 2M — h — 1 and summation over 
repeated indices is understood. The Lagrangian @ is invari- 
ant under the group 0(2TV). For N = 1, this reduces to the 
Lagrangian for the nonrelativistic Bose gas, and we can iden- 
tify the coupling constant with the s-wave scattering length, 
g — 8iva 2 . In order to eliminate the quartic interaction from 
the Lagrangian we introduce an auxiliary field a and add to 
Eq. © the following term 
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The Lagrangian can now be written as 
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By using the equation of motion for a, one can eliminate this 
field altogether and one recovers the original Lagrangian 

Integrating over the fields tj}2,ip3—ipN, we obtain the effec- 
tive action for a and tpv 



S cB = (N - l)Trln [d T - V 
+ / dr I d 3 X 
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where j3 = 1/T. In Eq. 0, Tr implies taking the trace of the 
differential operator inside the brackets. We next parametrize 
the fields ipi and a by writing them as sums of space-time 
independent expectation values 0o and im and quantum fluc- 
tuating fields tp and a: 

N(fo + $i , (6) 

= ■ (7) 

TV 



a — im + 



The rescaling with factors of 1/yN is just facilitate the count- 
ing of factors of TV. The field <f>o can be rotated so it is real. It 
can be shown that the expectation value of a is purely imag- 
inary so m is real. Substituting Eqs. into Eq. iJKJ, we 
obtain 



S cB = (TV-l)Trln 
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where we have set fcs = 1. Writing ipi = (<f>i + i<p2)/v2 and 
expanding Eq. (JHJ in powers of 1/vTV through order 1/TV, we 
obtain 
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where V is the volume of the system, \ = 
(4>i(P), 02 (-P), &(P)) are the Fourier transforms of cj>\, 
cj>2, and q. We have neglected linear terms that vanish at 
the minimum of the thermodynamic potential. We have here 
introduced the sum-integral 
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The integral over three-momentum is regularized by using a 
three-dimensional cutoff A. Finally, the function II(P, m) is 
defined by 



U(P,m) = jf 

Jc 



1 



q q 2 + 2m + iqo (p + q) 2 + 2m ■ 



i(po + qo) 

(11) 



1 In principle, one could introduce TV different chemical potentials 

but this is not natural in this context. The symmetry group 
would then be [0(2)]^. 

2 Note that the expansion parameter formally is 1/2N so the ex- 
pansion makes sense also for TV = 1. 



III. THERMODYNAMIC POTENTAL AND 
GAP EQUATIONS 

The thermodynamic potential Q is given by all one-particle 
irreducible diagrams and is given by a series in 1/TV. The 
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first two terms are easily obtained from Eq. by gaussian 
integration over \- We can then write 



A^LO + f^NLO , 
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where 
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NLO = 
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(14) 

The free energy .F is given by all connected diagrams and is 
independent of the condensate (j>o and the expectation value 
m. So if the thermodynamic potential is evaluated at the 
values of 0o and m, which satisfy the equations 

^ 0, (15) 
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the one-particle reducible diagrams vanish. The values of 0o 
and m, which satisfy the stationarity conditions 1151 and 116H 
can be written as series in 1/N: 



00 = 00, LO + 00, NLO 

1 

m = m L + -^TTlNLO ■ 



(17) 
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By series expanding the thermodynamic potential around the 
LO solution to the gap equations, the free energy can then be 
written as 

T = AT^LO^LO, <t>o, lo) + ^NLo(m.LO, 00, lo) 

d^LO 



+7TiNLO- 



+00, NLO 



dm 



d<f> 



(20) 

(21) 
(22) 



m=m L o ,4>0 = 4>a, LO 

+ C(l/iV 2 ) .(19) 

™=™lo ,<t>a=4>a, lo 
The free energy to NLO then reduces to 

T — NTlO + ^NLO , 

where 

.FlO = ^Lo(n"lLO, 00, lo) , 

■FnLO = ^NLo(«lLO, 00, Lo) • 

The number density p is given by the expectation value 

P = <VV) • (23) 
From the path-integral representation of the free energy T 

e- pVr = j Vll>*T>iPe- f ° dT 1 d3xC , (24) 

and Eq. the expression for the number density can be 
written as 

^ dfj, 



(25) 



Using Eqs. Q and (12511 . one finds the following exact result 
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A. Leading-order results 



At leading order the gap equations 1151 and 1161 1 become 



4m0o = 



, 2 p + 2m 

bo = 
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(28) 



The solution to Eq. I27H is that either m or <j>o vanishes, m = 
corresponds to the Bose-condensed phase and we discuss this 
first. 

Substituting m — into Eq. 1131 . the leading-order conti- 
bution to the free energy reduces to 

■^lo = + ^ In [ipo +P 2 ] ■ 

The sum-integral can be calculated by first writing the sum 
over Mastubara frequencies as a contour integral. By per- 
forming this contour integral and integrating over angles, we 
are left with 



2g ». J0 

2 A 5 



/'• 



5/2 



(29) 



The divergent term is eliminated by adding a vacuum coun- 
terterm AjT to the free energy T and the renormalized LO 
free energy is 
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Note that the LO free energy density does not agree with 
the standard mean-field result. The first term in Eq. (13011 is 
the zero-temperature result in the mean-field approximation, 
while the second term is the finite-temperature contribution 
for an ideal Bose gas. The LO number density in terms of the 
condensate 0o follows from Eqs. 1261 and I28K 
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The sum-integral in the gap equation I3H can be calculated 
in a similar manner as the one appearing in the free energy 
and one obtains 
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The divergent term reflects the fact that the expectation value 
of the number density operator also has divergences and we 
need to add a counterterm Sp to eliminate them [is?] . This 
yields 



P = 00 + 



3/2 
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(33) 



The critical temperature is found by setting 0o = 0. This 
yields 
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T" = 4tt 
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(34) 
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which is the standard result for an ideal Bose gas |15|. In- 
verting Eq. Ij45^ then gives the critical number density p° for 
Bose condensation 
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We next discuss the symmetric phase where <j>o = 0. The 
free energy in the symmetric phase is given by 



Hi 
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(36) 



where Li n (x) — J2 k x k /k n is the polylogarithmic function and 
we have have eliminated p, + 2m in Eq. (11311 using Eq. 1161 . 
In the symmetric phase, the number density reduces to 
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from 14011 and reads 
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Another more useful expression for the number density follows 
from Eqs. JUJ and 
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The 1/N expansion is not limited to a weakly interacting 
Bose gas. If we make the additional assumption of weak cou- 
pling, we can make further approximations by simply neglect- 
ing g!l(P, m). This approximation is valid outside the critical 
region and the free energy reduces to 



it 1 
2g T 2 



In [po + e 2 (p)] 



(43) 



where 



where we have neglected a temperature-independent diver- 
gent term. The number density in the symmetric phase can 
be interpreted as that of an ideal Bose gas with an effective 
chemical potential /i e ff = —2m. Eq. 1281 can be interpreted 
as a self-consistent one-loop equation for the chemical poten- 
tial and the critical number density for Bose condensation is 
then given by /Lt e ff = 0. 



B. Next-to- leading order results 

The NLO gap equations 1151 - 1161 are 

W I 
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where we have used that it is consistent to set m = in the 
NLO terms. Setting N — 1, the expression for the free energy 
now becomes 



T = - 



In [ip + p ] 

n(P,o ) + i + ^!_ 
g Po+p 



(40) 



where </>o satisfies the LO gap equation 1281 and is to be con- 
sidered a function of /i. The number density can be derived 



e(p) = p\Jp 2 + 2g4>l , 



(44) 



is the well-known Bogoliubov spectrum Similarly, the 
number density reduces to 
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After performing the sum over Matsubara frequencies, we ob- 
tain 



d a p e 2 (p)/p 2 +p 2 



(2tt)3 



[l + 2n(e(p))] , (46) 



where n(x) = l/(e^ x — 1) is the Bose- Einstein distribution 
function. At T — 0, this reduces to the Bogoliubov approx- 
imation and gives the standard weak-coupling result for the 
depletion of the condensate 0: 
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(47) 



At finite temperature, this approximation goes beyond the 
Bogoliubov approximation since we via Eq. 1451 take into 
account the temperature dependence of the condensate 4>o 
when evaluating 1431 [lr^|. 

In order to find the critical number density p c , we must 
evaluate the right-hand side of Eq. 1461 in the limit <f>o — * 0. 
The Bogoliubov spectrum then reduces to the free-particle 
spectrum and we recover the mean-field result 1351 . 

We have seen that by neglecting II(P, m), one recovers the 
mean-field results for the dilute Bose gas. In order to go be- 
yond mean field, we must keep the self-energy in the gap equa- 
tions. Since the phase transition is dominated by the nonper- 
turbative long-distance physics, the dominant contribution 
comes from the static Matsubara mode I3ll 1321 licl. The 



4 



contributions from the nonstatic Matsubara modes can be 
treated perturbatively and this fact underlies the calculations 
of Refs. [3l| EJJ using dimensional-reduction techniques. 
If we denote the contribution to the self-energy U(P,m) from 
the zeroth Matsubara frequency at po — by Ho(p, m), we 
have 



where we have used g — 8na. Thus the critical number den- 
sity decreases linearly with the scattering length a. Instead of 
calculating the critical number density at fixed temperature, 
we can calculate the critical temperature at fixed number den- 
sity by using 
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The critical number density is then found by settting 
and inserting 14911 into 1421 and rearranging terms, 
yields 
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(50) 



(51) 



can be interpreted as the leading correction to the self-energy 
function in the 1/N expansion. The corresponding string of 
Feynman diagrams is shown in Fig. The first of these dia- 
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which is valid to leading order in the number-density expan- 
sion [HIH. Thi s yields 
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It is interesting to note (2a that the results 152H and 154H 
are independent of N, although they are correct only in the 
limit N 00. The value 2.33 of the constant c compares 
reasonably well with c ~ 1.3 from lattice calculations [34. 135L 
13^ . an d classical field-theory simulations |3*3 | . 



IV. SUMMARY 




FIG. 1: Diagrams contributing to the self-energy at LO in the 
1/iV expansion. 

grams corresponds to the HF approximation. The self energy 
needs mass renormalization, which is carried out by subtract- 
ing S(0) |49| . Similarly, the second term in Eq. H50t can be 
interpreted in terms of Feynman diagrams 39, 40] . These are 
shown in Fig. [5] 




FIG. 2: Diagrams contributing to the number density at NLO 
in the 1/N expansion. 

The first term in 1501 is p°, while the the second term was 
first calculated by Baym, Blaizot, and Zinn- Justin [39|. and 
later by Arnold and Tomasik |40j . It is given by gT/32n 2 and 
the critical number density reduces to 



Pc = Pc 



T\' 2 
4- (p° c ) 1/3 g 



[<(§)] 



4/3 



(52) 



In this paper, we have used the 1/iV-expansion to derive 
thermodynamic quantities of the interacting Bose gas. If we 
make the additional assumption of weak coupling, we obtain 
the standard results for the dilute Bose gas. If we neglect 
the contribution to the number density from the nonstatic 
Matsubara modes at criticality, we obtain the result for the 
critical number density first obtained by Baym et al l3ll |3^| 
The 1/N expansion gives a coherent description of the dilute 
Bose gas in the critical region as well as outside it, where 
mean-field theories are normally sufficient. 

We have not calculated the critical exponents associated 
with the second-order phase transition of the Bose gas. How- 
ever, making the same approximations as we did for our calcu- 
lations of the critical number density, one expects to obtain 
the standard results for the 1PI effective action to NLO in 

i/iv m. 

It is natural to ask whether the 1/N expansion is reliable 
for N = 1. In the case of the critical temperature for Bose 
condensation, the 1/N correction was calculated by Arnold 
and Tomasik |4Cj| and the correction was approximately 26%. 
This is as expected: a correction of 1/2N or 50% multiplied 
by a constant of order one. 
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